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By introducing diffeomorphism and local Lorentz gauge invariant holonomy fields, we 

study in the recent article [S.-S. Xue, Phys. Rev. D82 (2010) 064039] the quantum Einstein- 
CD" . 

1 Cartan gravity in the framework of Regge calculus. On the basis of strong coupling ex- 

pansion, mean-field approximation and dynamical equations satisfied by holonomy fields, we 
•^T ■ present in this Letter calculations and discussions to show the phase structure of the quantum 

Einstein-Cartan gravity, (i) the order phase: long-range condensations of holonomy fields in 
strong gauge couplings; (ii) the disorder phase: short-range fluctuations of holonomy fields 
in weak gauge couplings. According to the competition of the activation energy of holonomy 
fields and their entropy, we give a simple estimate of the possible ultra-violet critical point 
and correlation length for the second-order phase transition from the order phase to disor- 
der one. At this critical point, we discuss whether the continuum field theory of quantum 
Einstein-Cartan gravity can be possibly approached when the macroscopic correlation length 
C*~) . of holonomy field condensations is much larger than the Planck length. 
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Introduction. Since the Regge calculus [l|, [2| was proposed for the discretization of gravity t 



he 

3 . 



5, 



ory in 1961, many progresses have been made in the approach of Quantum Regge calculus 
and its variant dynamical triangulations |(| . In particular, the renormalization group treatment is 
applied to discuss any possible scale dependence of gravity . Inspired by the success of lat- 
tice regularization of non-Abelian gauge theories, the gauge-theoretic formulation [8] of quantum 
gravity using connection variables on a flat hypercubic lattice of the space-time was studied in the 
Lagrangian formalism. The canonical quantization approaches to the Regge calculus in Hamilto- 
nian formulation are studied in Ref. [9]. All these studies are very important steps to understand 



the Einstein genera 
our recent articles 



relativity for gravitational fields in the framework of quantum field theory. In 



10, 



111 ], by introducing diffeomorphism and local Lorentz invariant (i.e., local 



gauge-invariant) holonomy fields, we present a diffeomorphism and local Lorentz invariant regular- 
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ization and quantization of Euclidean Einstein-Cartan (EC) gravity in the framework of quantum 
Regge calculus. 

Based on this theoretical formulation of quantum Einstein-Cartan gravity, in this Letter we 
present a preliminary study of the possible phase structure and ultra-violet critical point for the 
second-order phase transition of the theory. On the basis of strong couplin g ex pansion, mean-field 
approximation and dynamical equations satisfied by holonomy fields 



10, 



111 ], some calculations 



and discussions are presented to show the phase structure of the quantum Einstein-Cartan gravity, 
(i) the order phase: long-range condensations of holonomy fields in strong gauge couplings; (ii) the 
disorder phase: short-range fluctuations of holonomy fields in weak gauge couplings. Moreover, 
according to the competition of the activation energy of holonomy fields and their entropy, we give 
a simple estimate of the possible ultra-violet critical point and correlation length for the second- 
order phase transition. At this critical point, the minimal area (volume) element is shown to be 
the Planck one, in addition we discuss whether the sensible continuum field theory of quantum 
Einstein-Cartan gravity can be possibly approached when the macroscopic correlation length of 
holonomy field condensations is much larger than the Planck length. The possible relation of this 
macroscopic correlation length to the cosmological constant scale is also discussed. 
Simplicial manifold. The four-dimensional Euclidean manifold 1Z 4 is discretized as an ensemble 
of No space-time points (vertexes) "x € 1Z 4 " and Af± links (edges) "Z^(x)" connecting two neighbor- 
ing vertexes. The edge (1-simplex) denoted by (x,//), connecting two neighboring vertexes labeled 
by x and x + in the forward direction fi, can be represented as a four- vector field Z^(x), defined 
at the vertex "x" by its forward direction [i pointing from x to x + and its length 

a^(x) = \l(t(x)\ ^ 0, l^{x) = ae^(x), (1) 

where the fundamental tetrad field e^(x) = e^(x)7 a is assigned to each edge (1-simplex) of the 
simplicial complex, and a is a characteristic length of the simplicial manifold Ai{a). On the edge 
(x, fi), we place U^(x) = e l 9^^) ; an SO(A) group-valued spin-connection fields oj^{x) = uj c ^{x)a a } ) . 
The fundamental area operator of the anti-clock like 2-simplex (triangle) h(x) is defined as 

SU X ) = W*) S h (x) = \S%,(x)\. (2) 

The fundamental volume element around the vertex "x" is defined as 

dV(x) = dV h (x), dV h {x) = Sg(x) (3) 

h(x) 

where J2h(x) indicates the sum over all 2-simplexes h(x) that share the same vertex x. The 
characteristic length a is a running length scale, a\ — > 0,2 ■ ■ ■ ajy-i — > ajy and a\ > 0,2 ■ ■ ■ ajy-i > a^, 
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correspondingly simplicial manifold M.{a{) — > M.(a,2) ■ ■ ■ M.{cln-i) —> M(cln). In the sense of 
Wilson renormalization group invariance, we will try to find a physical scaling region where the 
macroscopic correlation length £ of the simplicial manifold is much larger than characteristic length 
a that is approaching the Planck length a^i = {SttG) 1 / 2 . 



Invariant holonomy fields. In Refs. 



10, 



111 ], introducing the vertex field v /J , u (e fl ,e l ,), we define 
the diffeomorphism and local gauge-invariant holonomy field along the loop C on the Euclidean 
manifold 1Z 4 



X c (v,uj) = Vctrexp <^ig j)^v fiu (x)u} fl (x)dx u j , (4) 

where Vc is the path-ordering and "tr" denotes the trace over spinor space. The regularization of 
the smallest holonomy field along the closed triangle path of the anti-clock like 2-simplex h(x), 

X h (v,U) =tr[v u ^(x)U fJ ,(x)v lip (x + a IJi )Up(x + a fl )v pu (x + a u )U u (x + a l/ )]. (5) 

The diffeomorphism and local gauge-invariant regularized Einstein-Cartan action, 

A EC = Ap{g,X h ) + A H (rt,Xh), (6) 

A P (e,U) = —2 [ x h(v,U)\ v ^ x ) = e lu ,( x y ys +h.c], 
g heM 

A H {e,Up) = —3- [ x h(v,U)\ Vfli/{x)=e ^ {x) +h.c] , 
y * heM 

where e fMU (x) = | [e p (x)e u (x) — e u (x)e fl (x)], the Immirzi parameter 7^0 Q, [3] and the gauge 
coupling g depend on the characteristic length a. In the naive continuum limit a — > a p \ and 
agujp <C 1: Eq. © approaches Einstein-Cartan action [see Section III(F) and Appendix B in 
Ref. [ll|]], when the running gauge coupling g(a) satisfies 

G e s = - g(a) — =^ G = — . (7) 
The partition function and the vacuum expectation value are defined as, 

Zec = J VeVUexp -A EC , <•••> = / VeVU{- ■ ■) exp -A EC - (8) 
The {Xh) obeys the dynamical equation, 

w = <M<f)>-<^J^f)>- w 

It should be mentioned that if O is not a diffeomorphism and local Lorentz gauge invariant operator, 
its vacuum expectation values must vanish (O) = 0, because diffeomorphism and local gauge 
symmetries are exactly preserved without any either explicit or spontaneous breaking. 
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Mean-field approximation. In order to show the phase structure and transition of regularized 
Einstein-Cartan theory ([6]), one needs to calculate (Xh) as a function of the gauge coupling g 
and the Immirzi parameter 7. The diffeomorphism and local Lorentz gauge invariant (Xh) acts 
as an order parameter. However, an analytical calculation of (Xh) is rather difficult for its non- 
perturbative nature. We adopt the mean-field approximation, though it is not diffeomorphism and 
local Lorentz gauge invariant, and try to gain some insight into the phase structure and transition 
of regularized Einstein-Cartan theory. 

Ill ], introducing the mean- field value ~ (v 2 ) and averaged area 



In Section VI of Ref. 

(Sh) = M^a 2 , for each 2-simplex h we define the local mean-field action Ah for the 2-simplex h(x) 



A 



tr 



ev(»)r^(a;)e /1 (a;) - e^(x)T^(x)e v (x) 



and the local mean-field partition function 



Z h = / VUVeexp-Ah, / VUVe = / dU^dU^dUpde^deu, 

J h J h J h 



where 



Ml 
8g 2 



75 
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U u (x)U p (x + a u )UUx) 



+ h.c, 



(10) 



(11) 



(12) 



Thus, the regularized EC action ([6]) and partition function (|8j) are approximated by their mean-field 
counterparts, 



heM heM 



(13) 



Eqs. (|10H13p are the mean-field approximation to the regularized Einstein-Cartan theory ()5|6l8p . 
In addition, we adopt the strong coupling expansion in powers of M^/Sg 2 to make analytical 
calculations. We approximately calculated the free-energy 



J-f c p (M h , 5 ,7,a) = -M\n(l + y h )-N 



2Vh 
1 + 2/h 



+ (Aec)o, Uh 



7 2 + l 
64c/ 4 7 2 d] 



\1. 



hi 



(14) 



where TV = J2heM ^ s ^ e total number of 2-simplexes and the mean-field value (Aec)o is an 
average with respect to Z^ c (|13p . By minimizing the free-energy (|14p with respect to Mh, we try 
to determine the mean-field value M£(g, 7), at which the free-energy (|14[) reaches its minimum. 



In Ref. 



Ill ], the local mean- field partition function and free-energy (|14|) have been obtained 



in the strong coupling expansion in terms of M^/8g 2 , up to the term (r^) 2 ~ 0[{M 2 y / Sg 2 ) 2 ] (see 



Eqs. (E.5) and (E.6) of Ref . [lj). 




A' 





B 



FIG. 1: We sketch a graphic representation of the dynamical equation ([9]) for the smallest holonomy field 
(Xh) (JSJ). Note that A and A 1 are the same vertex, so are B and B' . In the right-hand side of the 
graphic equation, the summation over all 2-simplexes h(l) associated to this edge is made. This figure is 
reproduced from Fig. 3 in Ref. [ill ]. 



In the previous article [ll|] , in order to estimate the minimal area of four-dimensional dynamical 
simplicial manifold, the term {Aec)o m the free-energy (|14p was very approximately calculated 
[see Eqs. (191), (192) and (193) of Ref. [11]]. In this Letter, in order to gain some sights into 
the phase structure of regularized Einstein-Cartan theory ([6]), we try to improve the calculation 
{Aec)o i n the free-energy (|14p by using strong coupling ex pan sion and the dynamical equation 
([9]). This is analogous to the mean- field approach developed TJ] for non-perturbative calculations 
of the Wilson loop in the lattice QCD. 

Dynamical equation. We try to solve the dynamical equation for the smallest holonomy fields 
in the framework of strong coupling expansion and mean-field approximation. The vacuum expec- 
tation value ([8]) can be written as 

(•••) = — ; — n r — r: — - = (• • -}o + hight - order terms, (15) 

(^-{Aec-Aec)) 

where hight-order terms stand for the series of strong coupling expansion (1/g 2 ) of exponential 
factor e~^ EC ~^ EC ^ in both nominator and denominator. Up to the leading order, replacing (• • •) 
by (• • -) in Eq. ([9]), we approximately write Eq. as follows, 

(X h ) « (X h ) ^Ag 2 (A EC )o/M 
(A E c)o « (A E c(u, 6 -^)) -(A EC (ul 6 -^)) , (16) 

where the first line bases on Eq. @. In below, we try to approximately calculate the right-handed 
side of Eq. (fT6"j) to obtain (Aec)o as a function of and 1/g 2 . Then, substituting the result 
(Aec)o = f(M^, 1/g 2 ) into Eq. (fT4"]). we obtain the approximate free-energy as a function of 
and 1/g 2 . 

The graphic representation of Eq. (fl~6j) is given in Fig. [TJ To obtain {Aec)o from Eq. ([TBI , we 
need to calculate the following four types of diagrams, shown in Figs. [2] and [3] . For the diagrams 
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represented in Fig. [21 using the mean-field approach (|10H13|) and indicating [• • •] to be anti-clock 
like and [• • -]t clock like, we have 



(FigEtleft)) « ]<tr[e v r^e M -e M r*ie v 



tr 



e v'^JfM> e Li' e V'^Jfi' e v' 



-\h,2 



1 



Zh 1 Zh 2 
x tr 



[det^/]/^ [dedU]h 2 exp - [A^ + Ah 2 ] 



tr 



~,h 2 



E (Pi ~ 

heM 



Zj ' 



(17) 



where in the last line we use Eqs. (E1)-(E7) in Appendix E of Ref. 111 ]. Due to the relations 
between anti-clock like orientation and clock like orientation r^L = r£t and e^r^e^ = — e^r^„e v 
(see Fig. [2]) , we have 



<Figj2£right)) 



•(tr j^e^r^e^ e^T^e^ 



tr 



>.} 



I o 



fc 2 



where in the last line we use Eq. (fT7|) . For these two cases, ^ //, v ^ v 1 and 



C ■/ - 



(2!±i) 2 + ±j(2i±i)- 2 , 
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(18) 



(19) 



and T>i is the total number of 2-simplexes h associating the fixed link I [19j. For the diagrams 
represented in Fig. [3l we have 



(FigE{left)} c 



— / [dedt/]^ exp- 



tr 



tr 



P r fol p — p r fel p 



4 



r , -i 2 

[dedU] hl tr(e u T Je^) exp-[Ai] 



3 E °r 

heM 



[dU^dU u dU p ] 

z h r 



-) 7 det [i-r^] 



(20) 
(21) 



and 



(FigEtright)) « (tr 



tr 



7 



[dedU] hl exp- [A hl ] 



x tr 



P r hi f= — p v hi p 



tr 



— J [dedU] hl [tr(e u r h ^)\ [tr^r^)] exp - [A hl ] 

- 2 S (t 



(22) 
(23) 



In lines flZQ and (g3J, we use Eqs. (183)-(185) and Eqs. (E5)-(E6) in Appendix E of Ref. [lit], and 
in the lines (f20|) and (f22j) . we use 

3 



J de^de v elel exp-A h = - r^jdet" 1 !/ - F f 



As a result, we obtain the (Aec)o as a function of and l/<? 2 , 

(A E c)o « AT |2(2?, - l)C r (j|) 2 + 3C r (I) + 2 



(24) 



(25) 




FIG. 2: Two simplexes /ii and /12 are not overlap and have a common 1-simplex (link) in the direction 
p. There are (T>i — 1) possibilities. Left: this is the second graphic representation in Fig. [TJ Right: this is 
the third graphic representation in Fig. [T] 





FIG. 3: Two simplexes h\ and hi are completely overlap. Left: this is the second graphic representation in 
Fig. [T] Right: this is the third graphic representation in Fig. [T] 

free-energy and two-phase structure. Substituting Eq. (|25j) into Eqs. (fTIj) , we obtain the 
approximate free-energy 



1 



JJ p c p (M h , g , 7 ,a) 



ln(l + y h ) + 2(V t - l)C r 



B- 

Zu 



+ 3CV 



(26) 
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Its minimum, as shown in Fig. 2] locates at 

M h * « yf n f (2 9 ) » 0.91 5 , (27) 

where y™ m ~ 0.04, the Immirzi parameter 7^1 and fundamental representation d-- = 4. The 
result ([27]) shows a nonvanishing mean-field value M£ decreases as the gauge coupling g decreases. 
The mean-field value (|27p shows that the framework of mean-field approximation and strong 
couping expansion is self-consistent for the expanding parameter (M^) 2 / (8g 2 ) ~ 0.1 being smaller 
than one. 

Submitting y™ m ~ 0.04, the location of free-energy minimum, into Eq. (|25p . we obtain the 
vacuum expectation value 

(A EC )/N « {A EC )o/M ~ 0.1. (28) 

From the first line of Eqs. Q, we approximately have 

(X h ) « (X h ) « 4 5 2 (^c) /AA ~ 0.4 5 2 . 

The mean-field values for the 2-simplex area ([2]) and the volume element ([3]) are 

(S h {x)} = a 2 M£ « 0.91 #a 2 , (dV(z)) = a 4 iV h (M^) 2 « 0.83 # 2 a 4 iV h , (29) 

where is the mean value of the number of 2-simplexes h(x) that share the same vertex. These 
nonvanishing values (]27H29p characterize an order phase in strong gauge couplings, as will be 
discussed below. 

Henceforth, Eqs. (|4I5|) will be called X-loop for short. The regularized Einstein-Cartan action 
Aec i s actually the ratio of the activation energy per area (the smallest X-loop Xh) and squared 
gauge coupling g 2 ("temperature"). In the order phase ((Xh) ^ or M£ ^ 0) for large coupling 
g 3> 1, the X-loops ([!]) are not suppressed because the smallest X-loops Xh undergo condensation 
by jointing together side by side to form surfaces whose boundaries appears as large X-loops. 
Namely, these X-loops proliferate and become macroscopic in the length £ that is the coherence 
correlation length of the system, leading to the area law 

(X C ) ~ ({X h )) n ~ exp[-A m i n (C)/<5 h (x))] J (30) 

where n = A m i n (C) / (Sh(x)) is the minimal number of 2-simplexes filling the minimal area A m \ n (C) 
that can be spanned by the loop C and (Sh(x)) ~ a 2 M£ is the averaged area of 2-simplexes. 
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0.15 




0.10 - 



1/6 



0.05 - 



4 



FIG. 4: In the strong coupling region (g > 1), the approximate free-energy (|26| is plotted as a function of 
Mh (a = 1). It shows that the location of minimal free-energy M£(g) becomes small, as the gauge coupling 
decreases, i.e., M^(gi) < M^{g 2 ) < M£(g 3 ) for gi — 1.4, g 2 = 1.5 and g 3 = 1.6. In the weak gauge coupling- 
region (g <C 1), a speculated free-energy with its minimum at M£ = is sketched. We approximately adopt 
T>i « 6 in Eq. 1)250 for a four-dimensional simplicial manifold, and results are not sensitive to T>i values. 

The result (|27p is obtained by solving the dynamical equation (|16|) in the framework of mean- 
field approximation and strong coupling expansion for (M£) 2 /(8g 2 ) < 1. Therefore, this result 
(|27j) does not apply to the weak gauge coupling region g <C 1, so that we cannot conclude — >• 
as p — >• 0. We have not so far been able to do any analytical calculation in the weak gauge 
coupling region of the regularized Einstein-Cartan action ([6]). Nevertheless, we can gain some 
insight into the possible phase in the weak coupling region by looking at the limit of gauge weak 
coupling g — > of the regularized Einstein-Cartan action ([6]). In the limit of weak gauge coupling 
g — >• 0, as we can see from Ap and Ah, the configurations of tetrad fields {e M (x)} and gauge fields 
{[/^(x)} have to be frozen to the configurations of small fluctuating fields for small Xh, otherwise 
the partition function (JSj) would vanish. Namely, tetrad fields {e^(x)} and gauge fields {U^(x)} 
undergo fluctuations at small scale a with large entropy. The smallest X-loops Xh are suppressed 
by their activation energy, and become then irrelevant for the large scale behavior of the system. 
Therefore, in the weak coupling region g <C 1, we conjecture the existence of the disorder phase 
with (Xh) = 0. In the framework of mean- field approximation, this means that the minimum of 
the free-energy locates at MjJ = 0, as sketched in Fig. HI 

These two distinct phases in the strong and weak coupling regions are characterized by the 
order parameter (Xh) or the mean-field value M£ in the framework of mean-field approximation. 
Taking into account the Immirzi parameter 7 > 1 in the action Ah [see Eq. (J6|)], we find from 
Eq. (|12p that the increasing value of Immirzi parameter 7 effectively leads the decreasing of gauge 
coupling g. Therefore we conjecture that the phase diagram should be the one sketched in Fig. [5l 
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0.06 
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0.01 
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(X h ) = 



M I ^ 
(Xh) o 

0.2 0.3 



1/5 



FIG. 5: We sketch the conjectured phase diagram in terms of inverse gauge coupling 1/g and Immirzi 
parameter 1/7. As discussed in the text, the critical point indicated g c — 4/3 for 1/7 = 0. While the critical 
line and point 1/g = and I/7 — 0.068 are arbitrarily sketched for indicating two-phases structure. 

Phase transition and critical coupling. Since we have not so far been able to calculate the 
order parameter (Xh) in the weak coupling region, we cannot exactly determine the critical point 
or line of the second-order phase transition from the order phase to disorder phase. Nevertheless, 
in this section, we try to discuss the critical point or line of the second-order phase transition. 

As indicated in Fig. 0J for the disorder phase in the weak coupling region, the minimal free- 
energy is zero locating at (Xh) = 0; for the order phase in the strong coupling region, the minimal 
free-energy is negative locating at (Xh) 7^ 0. The second-order phase transition from the order 
phase to disorder phase occurs when the minimum of free-energy for weak couplings flips into the 
one for strong couplings, as a result of the competition between the activation energy of X-loops 
and their entropy. We expect the second-order phase transition taking place at a critical coupling 
g = g c 7^ for 7>1, and we try to estimate it in a simple way that was adopted to estimate the 



critical 
helium 



joint of the second-order phase transition for two-dimensional systems [ijj, the superfluid 



161 ] and the U(l) lattice gauge theory [1 



For this purpose, we consider the partition function of a single AT-loop of arbitrary length C 
given by the integral 



Z c = J VeVU(X c )exp-A E c 



exp 



1 



(25) 2 



(31) 



heA(C) 

where (i) JT>A(C) is the functional measure of all possible surface-area A(C) bound by the closed 
loop C; (ii) (2d) A ^^ Sh ^ approximately accounts for the number of possible configurations (surface 
deformations at the area scale (<Sh)) of the surface with a given surface-area A(C), this is related to 
the entropy of the surface-area A(C); (iii) (Xh) is the activation energy of a 2-simplex surface-area 
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(Sh), J2heA(C) (Xh) stands for the activation energy of the surface-area A(C) and (2g) 2 plays a role 
of "temperature". The number of local deformations of a 2-simplex area in a three-dimensional 
simplicial manifold is 6, and the number of local deformations of a 1-simplex length in a two- 
dimensional simplicial manifold 4. We assume this number to be 2d in a d-dimensional simplicial 
manifold. This assumption is not crucial, as you will see below, for a simple estimate of the critical 
coupling. The free-energy of a grand-canonical ensemble of arbitrary surface-area A(C) of the loop 
C is then given by Fq ~ — In Zq , and 



Z C ~ f VA{C) exp + -H \n(2d) -—2 £ W 
J (5 h ) (2 5 )2 ^ 

A(C) 



y i>i4(c) 



exp +- 



This integral converges only below a critical gauge coupling g, 



1 

(2s) 5 



ln(2d) - j^iXh) 



(X h 



1/2 



(32) 



(33) 



.ln(2d). 

At the critical coupling g c , a single 2-simplex configuration is activated and its activation 

energy should be the order of the Planck scale l/a p \, i.e., (X^) ~ 0(1), as the characteristic length 
of simplicial manifold is approaching the Planck length (a — > a p i). Eq. (|33p gives the critical 
coupling g c ~ O(l). 

Above the critical coupling g c , the integral diverges and the ensemble undergoes the second- 
order phase transition in which the surface-area proliferates and becomes macroscopically large 
with the coherent correlation length £ 3> a. This means that the surface-area A(C) of an X-\oop 
(Xc) can only be easily deformed beyond the scale £ 2 , indicating the "condensation of X -loops ", 
and the ensemble stays in the order phase. The divergent integral (|32[) can be written as 



Z c ~ J VA(C) exp + A(C)/f, 



for g > g c . Below the critical coupling g c , none of 2-simplex configurations (X^) is activated and 
{Xh) = 0, there are short distance (~ a) fluctuations of tetrad fields e^(x) and group- valued fields 
Ufj,(x) with large entropy, and the ensemble stays in the disorder phase. This implies that the 
second-order phase transition from the order phase to disorder phase takes place at the critical 
gauge coupling g c . 

In the order phase, as the characteristic length of simplicial manifold is getting smaller and 
approaching the Planck length (a — > a p i), the gauge coupling g(a) is approaching to the critical 
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gauge coupling g c , which is an ultra-violet fix point. The coherent correlation length £ (|34p becomes 
macroscopically large 

(S h ) i 1/2 



21n(2d) 
0.48 



1/2 



(9 -9c) 1 / 2 



> a, 



,-, ~ . (35) 
[9 -9c) L/2 

in the neighborhood of the critical coupling g ~ g c + + , where a quantum field theory of the 
Euclidean Einstein-Cartan gravity can possibly be realized. In the second line of Eq. (|35p . d = 4 
and the mean-field value of 2-simplex area ([29]) are used. In Eq. ([35]) . the critical coupling g c ~ 0(1), 
critical exponent v = 1/2 and proportional coefficient Co ~ 0.48 are preliminary results obtained 
in this simple estimation. 

In the mean-field approximation we have (Xh) ~ (M^) 2 . (Xh) ~ 0(1) implies M£ ~ 0(1) in 
the nontrivial continuum limit (|35j) for a — > a p i and g — > g c . Therefore, the minimal averaged area 
of 2-simplexes in the mean- field approximation is given by Eqs. (|27p and (|29p 

a 2 M^ = 0.91#(5)5 2 



<S h (a)) 



a->-o p i,g-*-3 

which shows the minimal area (volume) element of the space-time is the order of the P 



(36) 
anck scale 



in the nontrivial continuum limit (|35p , the basic arena of physical reality we live on 

In addition to the case in the nontrivial continuum limit (|35p . we try to discuss this critical 
coupling g c by looking at the naive continuum limit of the regularized Einstein-Cartan action ([U]), 
where the gauge coupling g = g(a) depends on the characteristic spacing of simplicial manifold a. 
In the naive continuum limit a — > a p \ and gduj^ <C 1, we obtain the effective Newton constant [see 
Eq. 0], 

G B ff(S) = ^(a)^ ) (37) 

which has to approach to the Newton constant G = a 2 j/(87r) in the continuum Einstein-Cartan 
theory. This leads to g(d)\a^a pl — > 4/3 + + and the critical coupling g c = 4/3. 

We turn now to a general discussion. As the running gauge coupling g{d) is approaching to 
its ultra-violet critical point g c (g — > g c ) for a — > a p \, physical and dimensionful quantities m(g,d) 
should satisfy the renormalization group invariant equation, 



dm 



dm 



(3(g) 



dm 



0, /%) 



- dg(d) 
-a- 



(38) 



da da dg ' da ' 

because the coherent correlation length £, the physical scale, becomes much larger than a (£ S> a). 
The running gauge coupling <?(a/£) > g c can be expanded as a series , 

g(d/0 = g c [l + a (a/0 1/v + O[(a/0 2/l/ ]\ , (39) 
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leading to the /3-function 



(5{g) = -5 



~dg(a/0 



Po + Pi(g - g c ) + 0[(g - g c ) 2 ] 



(40) 



da 



where (3q = and /3i = —If v. Assuming m = £ , as the solution to the renormalization group 
invariant equation (|38p . one obtains that the coherent correlation length £ follows the scaling law 



where the proportional coefficient Co = (aoffc)^ and critical exponent za Eq. ([41 p has the same form 
as Eq. (|35p , Non-perturbative calculations by numerical simulations are required to determine the 
proportional coefficient cq and the critical exponent v in Eq. (|35[) or (|41[) . 

Some remarks. In this Letter, we present an analytical study of phase structure and critical 
point of the quantum Euclidean Einstein-Cartan gravity. For the order phase, calculations and dis- 
cussions are based on the approaches of strong coupling expansions, the mean-field approximation, 
and the dynamical equations for holonomy fields. For the disorder phase, we have not been able 
so far to do analytical calculations in weak gauge couplings region, the discussions on this phase 
are based on the limit case of gauge coupling g — > 0. The possible ultra-violet critical point and 
correlation length for the second-order phase transition are estimated in a simple model, according 
to the competition of the activation energy of holonomy fields and their entropy. Therefore, these 
results and discussions on the order- and disorder-phase structure, the ultra-violet critical point and 
correlation length for the second-order phase transition are preliminary. Numerical simulations are 
essentially required to check these preliminary results on the phase structure, ultra-violet critical 
point and correlation length for the second-order phase transition before one can conclude that a 
sensible continuum field theory of the quantum Euclidean Einstein-Cartan gravity can be defined. 
The coherent correlation length £ (|35p or (j41|) is an intrinsic scale of the quantum gravity, analo- 
gously to the intrinsic scale Aqcd of the quantum chromodynamics SU (3) theory (QCD) for the 
strong interaction. We speculate that the scale l/£ 2 might have some relation to the cosmological 
constant Acos> an d we leave this topic to a further work. 
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